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Abstract. If the polaron coupling constant a is large enough, bipolarons or 
multi-polarons will form. When passing through the critical ct c from above, 
does the radius of the system simply get arbitrarily large or does it reach a 
maximum and then explodes? We prove that it is always the latter. We also 
prove the analogous statement for the Pekar-Tomasevich (PT) approximation 
to the energy, in which case there is a solution to the PT equation at ct c . 
Similarly, we show that the same phenomenon occurs for atoms, e.g., helium, 
at the critical value of the nuclear charge. Our proofs rely only on energy 
estimates, not on a detailed analysis of the Schrodinger equation, and are very 
general. They use the fact that the Coulomb repulsion decays like 1/r, while 
'uncertainty principle' localization energies decay more rapidly, as 1/r 2 . 



1. Introduction 

We investigate the binding-unbinding transition of the ground state of multi- 
polaron systems, both in the original Frohlich model 6 and in the Pekar-Tomasevich 
large coupling approximation [18] . The parameter that varies is the strength of the 
electron-electron Coulomb repulsion, conventionally denoted by U > 0. (We could, 
alternatively, let the field coupling constant a vary, but it is mathematically con- 
venient and equivalent to vary U .) 

The conclusion of our rigorous analysis is that this transition is always 'first- 
order' in accordance with a proposal of Verbist, Peeters and Devreese [22]. That 
is, the Coulomb repulsion jumps discontinuously from a positive value to zero and 
the radius, too, jumps discontinuously. 

The same thing happens for the helium atom unbinding as U varies, as shown 
by M. and T. Hoffmann-Ostenhof and Simon [5]. Our analysis does not proceed, 
as in by analyzing the Schrodinger equation, but, instead, is based on energy 
estimates. Indeed, our proof is also valid in the presence of a magnetic field; the 
proof in [9] cannot accommodate a magnetic field because it relies on the positivity 
of the ground state wave function. To illustrate our approach we first show how 
to reproduce the result of [9] as a warm-up to the harder polaron problem. The 
overall lesson is that this kind of discontinuous binding will occur whenever the 
net repulsion at large distances falls of slower than r~ 2 . This conclusion does not 
depend on Bose or Fermi statistics. 
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Frohlich's Hamiltonian for N particles interacting with the longitudinal distor- 
tions of a polar crystal is 



N 

H u N) = E (p? - v^o*)) + ff / + (i.i) 



with X = (x l5 . . . ,xjv) G K 3W and 



^W=Vr--^- (1-2) 



i<j 



The Hilbert space is L (R d ) ® J 7 , where is the bosonic Fock space for the 
longitudinal optical modes of the crystal, with scalar creation and annihilation 
operators a'(k) and a(k) satisfying [a(k),a'(k')] = S(k — k'). The phonon field 
energy is 

H f = I dka\k)a(k) (1.3) 

and the interaction of the crystal modes with the electron is 

i r fjh 

4>{x) =-— — (e ik *a(k) + e~ ik ^(k)) , (1.4) 
V2vr J R 3 \k\ 

with coupling constant a > 0. We define the ground state energy to be 

E^ia) = inf specffW = ^ inf ^IH^I*) (1.5) 

and the break-up energy to be 



(7 



(a)= min ( (a) + E { V N ~ n) (a)) . (1.6) 

It is always the case that E^\a) < E^\a). If there is strict inequality, we say 
that the iV-particle system is bound. 

A note on particle statistics. As we said, the results in this paper do not depend 
on Bose statistics (restrictions to symmetric spatial functions) or Fermi statistics 
(restrictions to anti-symmetric functions of space and spin). The numerical values 
of the binding energies and the discontinuities, etc., will depend on statistics, but 
not the qualitative features proved here. 

Let us start by recalling some previous results. Some of these results are only 
valid in the Pekar-Tomasevich (PT) approximation, which will be described later. 

(1) Physically, U > 2a, but the case U < 2a has been considered in [8]. 
They show that for fermions the energy satisfies —C\N 7 / 3 < E^\a) < 
~C 2 N 7 / 3 . The same arguments for bosons lead to ~C[N 3 < E$f\a) < 
— C' 2 N 3 . In [I] [2] the sharp asymptotic constants for bosons and fermions 
are identified in the PT approximation. 

(2) The case U = 2a is special. In the context of the PT model [5] showed that 
the energy behaves linearly in N for fermions and in [5] an 7V 7//5 behavior 
for bosons is proved. 

(3) For U > 2a the thermodynamic limit, lini/v^oo Ey \a)/N, exists [S]. More 
importantly, it was shown in [5] that there is a value U*(a) such that 

J u 

Uli-l/U' T 111 f ! W i ll:-|1 .ll I/Ill 1 



E [ v\a) = NE {1 \a) for all U > U*(a) and all N. (For N = 1, we drop 
the index U in the notation E^\a).) In other words, there is no binding 



of any kind for U > U* (a) 
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(4) There is, in fact, binding for some range of a, U and N with U > 2a. (If 
N = 2 and a is large then binding occurs for 2a < U < 2.3 a [2T | fl9 | [23].) 

For every fixed N and a there will generally be certain critical values of U, which 
we generically denote by U c , at which binding will either appear or disappear. That 
is, 

E$\a)=E%\a) 

and 

El N \a)<E^(a) 

for either all U € (U c — S, U c ) or all U £ (U c , U c + S), for some S > 0. The usual case 
is that there is a unique U c such that binding occurs for U < U c and no binding 
occurs for U > U Cl but we do not know if this is true and we certainly do not know 
how to prove it. The problem is that the energies of the break-up components also 
depend on U, so there is no obvious monotonicity. 

The question we address here is what happens at U c , or as U approaches U c from 
the binding side. Does the compound simply explode into two pieces like a super- 
nova or does it get larger and larger like a red giant? Our answer is that it always 
explodes, that is, the maximal distance between the particles stays finite. For the 
bipolaron this was seen variationally in |22j and stated as a general mathematical 
question in |16j . 

One would like to say that this implies the existence of an eigenfunction when 
U = U c , but the polaron problem presumably never has bound states because of 
translation invariance [7]. (No doubt there are bound states of total momentum 
equal to zero, but we will not explore this here.) The PT problem, on the other 
hand, does have honest energy minimizers in a one-sided neighborhood of U c 
111)] , and we use the result about the finite maximal distance to show, in a separate 
theorem, that the PT problem has a minimizer when U — U c . 

In any case, a first order phase transition occurs. Suppose there is a break-up 
into two complexes of n and N—n particles each. Suppose particle one goes into the 
first complex and particle two goes into the second. Then their mutual Coulomb 
repulsion is zero after the break-up, but it is uniformly (in U) bounded away from 
zero before the break-up. 

While in a physical situation one is usually not at the critical value of the param- 
eters, the results proved here can be useful for a general understanding of quantum 
mechanics and, perhaps, for numerical analysis close to the critical values. 



2. The helium problem 

As a warm-up, we consider the ground state of a two-electron atom, of which 
helium and the hydrogen ion are examples. The Hamiltonian 

2—1 1 1 

acts on L 2 (M 6 ). The ground state energy E\j of Hjj is monotone increasing with 
respect to U, and there is a critical U c > such that for U < U c , Ejj is lower than 

-- = inf spec (p 2 - Ix]- 1 ) , 
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which is the energy of a hydrogen atom, and for U > U c they are equal. Our general 
method, which will later be applied to the polaron problem, is illustrated by the 
following theorem of M. and T. Hoffmann-Ostenhof and Simon [3]. 

Theorem 1 (Binding for helium j9j). Hjj c has a ground state eigenf unction ifju c 
in L 2 (R e ), i.e., H Uc ^ Uc = -\ij} Uc . 

Remark 1. The proof uses only two facts about the operators p 2 and | or | 1 : First, 
p 2 is non-negative and, second, the localization errors for p 2 fall-off faster than 
| a; | 1 . The proof would work for many other combinations having this essential 
relationship. For instance, as mentioned in the introduction, p 2 could be replaced 
by (p + A(x)) 2 for a suitable magnetic vector potential A(x). 

Theorem [T] follows by standard weak convergence arguments from the following 
Proposition [TJ The details of the arguments, showing that the ff 1 (R 6 )-weak limit 
of ground states ipu for U — > U c is not zero, are worked out in Section [5] on the PT 
equation, which is more complicated, in fact. Once one knows that the weak limit 
is not zero it is easy to see that this weak limit has to be a ground state. 

We know from [3] that U c > 1, and therefore for the proof of Theorem [1] the 
assumption U > 1 + 8 in the following proposition is not really a restriction. The 
main point is the uniformity as U — > U c . 

Proposition 1 (Upper bound on the helium radius). For any 6 > 0, there is 
a constant C'$ > such that for all 1 + 8 < U < U c and for all normalized ground 
states tpu of Hjj one has 

WuWx&Wu) >C Sl (2.1) 

where \x\oo = max{|xi|, |a; 2 |}. 

Proof. Our goal will be to prove the following operator inequality. (Actually, the 
analogous statement for expectation values in the ground state would suffice for 
our purposes.) There is a radius £o > and constants c, C > such that for all 
larger radii I > Iq and for all U > 1 + 8 one has 

Hu-E v > - J 6(t - |ar|oo) + (~ - E v + flfla^ - i) . (2.2) 



Here, 9 is the Heaviside function, i.e., 9(t) = 1 if t > and 9(t) = if t < 0. 
Since Ejj < —j, this inequality, evaluated in a ground state ipu, implies 

^ / ipfr dx > c / dx for all £ > £ a . 

Lemma [TJ stated below, implies (I2.ip . and our proof is complete. 

The crucial inequality (|2.2[) has a simple physical interpretation. The first term 
is a slightly negative energy, but vanishes relatively rapidly like £~ 2 . This leaves us 
with the second term, which is like a potential barrier. Although it goes to zero as 
£ — > oo, it does so much slower than £~ 2 'on average'. The significance of this is 
quantified by the following lemma. 

Lemma 1 (Calculus lemma simple version). Let p £ L 1 (R+) be non- 
negative, with f Q p(r)dr — 1. Assume that there are constants b > and £ c > 
such that 

b r e r°° 1 

p J o p(r)dr > ^ -p(r)dr (2.3) 
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for all £ > 4 . Then 



l p{r)dr * WTQ {2A) 



P(r)dr>j^. (2.5) 



A generalization of this lemma is given later in Lemma [4] 
Proof of Lemma [H We start by showing that 



-p(r)dr > — L- - I f C p(r)dr . (2.6) 

To prove this, we define f(R) = p(r) dr and rewrite the assumption as b£~ 2 f(£) > 
r^ 1 p(r)dr for all £ > £ c . Therefore 

1 f°° 1 f°° 1 fC/ 1 

ip(r)dr = / -f'(r)dr = / -f(r)dr - 

r Jo r Jb r z £ r . 



OO rOO 



>6- x / / Ip( s )dsdr-^ 



6- 1 (l-f(£ c )-£ c J™ Ip(a)dsJ 



/(4) 
4 



This is the same as (|2.6[) . 

In order to derive (|2.4I) from (|2.6|) we distinguish two cases according to whether 

C 1 Jo 4 P( r ) dr < (2(6 + 4))" 1 or not. In the first case, (gU) yields 

1 f°° 1 1 

—p(r)dr > / —p(r)dr > 



2(6 + 4 



In the opposite case, we have 



i i i r c i 

—p(r)dr > / —p(r)dr > — / p(r)dr > 



4 Jo " 2(6 + 4) ' 

This proves 

Finally, to prove (I2.5[) we combine the assumption (|2.3p at £ = 4 with (12.61) to 
get 



b r tc , , , r°° i , , . i i 



^ (r)dr ^X c r p( - r)dr ^bTT c -£- c J p{r)dr - 



which is the same as (12.5[) . □ 
We now turn to the proof of the potential-theoretic inequality (I2.2[) . 

Lemma 2 (Partition of unity). For any < e < 1/2 i/iere is a C e > swc/i i/iai 
f/ie following holds for any £ > 0. There is a quadratic partition of unity 

Xo(x) 2 + Xi(x) 2 + X2(x) 2 + X3(x) 2 = 1 for all x = (x u x 2 ) e x M 3 , 
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such that 

Xo( x ) = unless \x\oo < t, 

Xi{x) = unless \x\oo > 1/2 , \x\oo < (1 — e)|afi - x 2 \ , 

\x\oo>i/2, \x\oo > (l-2e)|a;i - x 2 \ 

\x\oo>£/2, |zU>(l-2e)|xi-z 2 | 



X"2( x ) = unless 
Xz{x) = unless 



\xi\ < (1 + e)|ar 2 ] , 
\x 2 \ < (l + e)N, 



and 



3 



< 



9i 



< 



c. 



l\x[ 



when Xo(x) 2 > , 

when xi(x) 2 + X2(x) 2 + X3(x) 2 > . 



One can choose xo an d Xi to be symmetric functions of (xi,x 2 ) and one can 
choose X2( x 1t x 2 ) — Xsi^^i)- This preserves both Bose and Fermi statistics. 

Sketch of proof of Lemma\^ The localization error . |Vxj | 2 is supported in three 
regions, {1/2 < ^ < £}, {(1 - 2e)\ Xl - x 2 \ < \x\ x < (1 - e)\ Xl - x 2 \} and 
{(1 + e) _1 |x 2 | < < (1 + e)|x 2 |}. It is a geometric question to check the stated 
inequalities on the supports of the different x's. We leave this to the reader. The 
fact that C e is independent of i is a consequence of scaling. □ 



Proof of inequality (|2.2p . By Lemma [5] and the IMS localization formula we can 
write, for any wave function "J, 



(*|H P |*)=2(*3 

3=0 \ 



fe=0 



X>ll^l 

3=0 



with ^j{x\,x 2 ) — ^(xi, x 2 )xj( x i, x 2 ) and with numbers Cj (depending on ^j). In 
the following, we shall derive lower bounds on ej . 

For j = 0, we simply use the bound on the localization error from Lemma [5] to 
bound 

C 

e >Eu--^. 

For j = 1, we use the fact that on the support of xii if, say, \x 2 \ = \x\oo then 

e e 



si > 



1 - e 



x 2 



1-e' 



Therefore, for any A > 0, 
-Nil" 1 - \x 2 

and 

ex > -2(e" 1 
For j — 2 we bound p\ — \x\\ 

so that 



> -e-'lx^ 1 > -2(e" 1 + \)lr l + X\x\^ 
A - CJ- 1 



-2 



*1 



> — j. Moreover, in that region 
and |xi - x 2 \ < (1 - 2e)^ 1 | 



e 2 > — 



1*2 



* 2 



U(l - 2e) - 1 - e - 



-l 
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Similarly, for j = 3, 



e 3 > 



1*3 



*a 



t/(l - 2e) - 1 



* 



In order to specify the free parameters e, £ and A we restrict ourselves to U > 1+5 
with some fixed 5 > 0. (We will not reflect the dependence on 5 in our notation.) 
First, we choose e so small that {7(1 — 2e) — 1 — e > 5/2 and we choose A = 5/2. 
Then we choose £ so large that 2(e _1 + A)^ 1 < 1/4 and C^q 1 < 5/4. With these 
choices, we have 



e, > 



5/4 



*, 



for any £ > £o and any j — 1,2,3. This, together with the fact that 

X0(xf < 0(£ ~ Woo) 

implies ([2~2]) . 

This completes our proof of Proposition [TJ 



□ 
□ 



Remark 2. One can ask a similar question for TV electrons and whether there is 
binding at the (unique) U — U c such that 



(2.7) 



(Here is the ground state energy of n electrons with repulsion strength U .) 
We denote by M the smallest number n < N — 1 such that E^ — E^ and 
we recall that by Zhislin's theorem [24] we have U c > M _1 . Whether equality or 
strict inequality holds has to be decided by an independent variational calculation 
(as Bethe [3] did for helium). We have nothing to say about the case U c — M _1 . 
On the other hand, if U c > M^ 1 our method should work and show that there 
is a bound state for the TV-electron system. Again, the method should extend 
to magnetic fields, pseudo-relativistic models with positive mass, spin-polarized 
systems, etc. 

Another area in which to try to utilize our method is Hartree or Hartree-Fock 
theory, though we have not pursued this. 



3. The Bipolaron 

We now return to our main theme and consider the simplest interesting case, 
namely two polarons, whose Hamiltonian is 

H { u ] = P \+pI- vW(zi) ~ ^<P{X2) + H f + U . (3.1) 

We have shown in [S] that there is a critical U c (a) < oo such that Efj(a) = 
2E^(a) for all U > U c (a), while E { J } (a) < 2E^(a) for all U < U c (a) (by 
concavity). It is easy to see that U c (a) > 2a. In the following, we assume that 
U c (a) > 2a. This is true for large a, since lim inf a,-^ U c (a)/a > 2.3 due to 
the convergence to the Pekar-Tomasevich functional in the strong coupling limit 
[SOUIIS], for which the critical ratio is known to be at least sa 2.3 [21 ] \19 \ [23 ] . 
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Theorem 2 (Upper bound on the bipolaron radius). For any e > 0, there is 
a constant C t > such that for all < 2a(l + e) < U < U c (a) and for all states \f 
one has 



1 



\xi - x 2 \ 



> 



U - 2a(l + e) 



2E^{a) - H 



(2) 



C e {l + U/a) 2EW(a)-E ( J\a) 

We emphasize two consequences of Theorem [5J 
(1) In any approximate ground state, in the sense that 



(3.2) 



^li^l*) < {l-\)2E {1 \a) + \E ( u\a) 



(3.3) 



for some A > 0, the expectation value of \x\ — x 2 \~ 1 is uniformly bounded 
from below by a positive number as U increases to U c {a). In particular, 
the size of the bipolaron does not increase indefinitely as U — ?> U c (a). 

( 2) 

(2) The bipolaron energy Efj (a) is not differ entiable in U at U — U c (a). 
While the right derivative is zero, the left derivative is at least as big as 
(U c (a) - 2a(l + e))/[C e (l + U B (a)/a)], for any e > 0. 

Proof. The proof of Theorem [2] will be divided into three parts. 

Step 1. Partition of the interparticle distance. As in [5] we choose a quadratic 
partition of unity and localize the particles according to their relative distance. In 
order to construct this partition, we pick some parameters b > 1 and £ > 0, and let 



<p(t) 




for t<l/b, 
™\¥m iov£/b<t<£, 
ioii<t<bt, 
for t>M. 



2 

TT t-l 

2 bl-l 



(3.4) 



For j > 1, let tpj(t) := ^(b^H), and for j = 0, let 



for t<£/b, 
for £/b<t<£, 



(3.5) 



for t > £ . 



Then 



^(pj(ty = l foralH>0. 

j>0 



(3.6) 



Using the IMS localization formula, we can write, for any wave function ^P, 



<*ii#>i*>=x;(^ 

3>0 \ 



hP-2J2H(\^-^\ 

k>0 



/ j>o 



with i&j(xi,X2) — ^(x\, X2)ifj(\xi — x 2 \) and with numbers ej (depending on 4^-). 
In the following, we shall derive lower bounds on ej. For our bounds we shall use 
the fact that on the support of <fij(\xi — x 2 \), the localization error is dominated by 



k>0 



< 



4(£-£/b) 2 



1 ifj = 0, 

&2(l-i) if j > 1 . 



(3.8) 
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Step 2. Estimate of localized energies. Because of p.8j) we can simply bound 



(3.9) 



J(J 2(£-e/b) 2 

For j > 1, we further localize each of the two particles to its own ball of radius 
b° L for some parameter L > 0. This will entail an additional localization error. 
Concretely, let 



1 



'27r|x| 



sin(7r|a?|) for |x| < 1 , 
for \x > 1 , 



(3.10) 

} 

and note that / dx\(x) 2 = 1 and J dx\Vx(x)\ 2 — k 2 . With 

^^u 1 ,u 2 (xi,x 2 ) = ^j{xx,x 2 ){VL)- i x{^ i {x l -u 1 )/L)x{b- j {x 2 - u 2 )/L) (3.11) 
we have, by a continuous version of the IMS localization formula, 

2tt 2 



H, 



(2) 



Vj) = / dui / du 2 (^j,u u 

1 JR 3 JR 3 \ 



H, 



(2) 



b 2 3L 2 



(3.12) 



Note that since VI > \xi — x 2 \ > W 2 £ on the support of ipj, the wave function 
j,«i,«2 is non-zero only if the distance d between the two balls of radius V L 
centered at u\ and u 2 , respectively, satisfies 



deDj := [b j - 2 e-4b j L,£b j ]. 



(3.13) 



We shall choose L < £/(4b 2 ). 



Lemma 3. Assume that \1/ is normalized and supported in B\ x B 2 where B\ and 
B 2 are disjoint balls of some radius R, separated a distance d. Then 

2a U 

1 

d 



(*l4 a >|*)>2EW(a)-^ 



(3.14) 



d + iR 

Eq. (|3.14p is an easy consequence of the functional integral representation of the 
ground state energy. The proof can be found in [SJ Lemma 1]. 

We shall apply inequality (|3.14[) to (13.121) . with U replaced by 2a(l + e) for some 
e > 0. Using also (|3.8j) we conclude that 

£/-2a(l + e) 



ej > 2E {1) (c 



2a 2a(l + e) 
d^Dj I. d + d + AWL 



\x\ - x 2 

- b- 2 > 



6 2 7T 2 



2{£-£/b) 2 



27T 2 

L 2 



(3.15) 



For given e > (and given b > 1), we can choose L to be a small enough constant 
times £ such that, as long as £ is large enough, the sum of the terms in the second 
line of (|3.15[) is positive. More precisely, for given e > (and given b > 1) we can 
choose 6 > small enough such that 

d + d + AS 



mm 

b- 2 -i6<d<l 



> 0. 



(3.16) 



With L = S£, we see that 

2a 2a(l 



for all j > 1 if 



d + AVL 



£>£c-= 



- b- 2 > 



2^2 



b 2 TT 



2(£~£/b) 2 



27T 2 

u 



> o 



4(1 -b- 1 ) 2 



1 

M 2 



(3.17) 
(3.18) 
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Under this condition, we thus have 

ej >2^( 1 )(a) + ||^|r 2 < 

for all j > 1. 



U-2a(l + e) 



\xi - x 2 \ 



(3.19) 



Step 3. Upper bound on the particle distance. From (|3.7[) and our bounds (I3.9P 
and (|3.19l) on ej we conclude that, for all £ > £ c , 

^2 



(2) 



2E ( - 1 \a) > fjS^a) - 2E (1 \a) 
U - 2a(l + e 



\xi - x 2 

\a)-2 

U - 2a(l + e) 
\xi - x 2 \ 



2£ 2 (l-b- 1 ) 2 
(1 - ip (\xi - x 2 \) 2 ) 



<po(\xi - x 2 \y 



>{E%\a)-2E^\a) 



2£2(1_ 6-1)2 

6{\x x -x 2 \-l). 



\xi - x 2 \ 



Now let "J satisfy 

(*|iJ^ 2) - 2E (1 \a)\q) < A (^ 2) (a) - 2E (1 \a 
for some < A < 1. For such a VP, the previous inequality yields the bound 

($> (l-\)(2E^(a)~E { J\a)^ + 



(3.20) 
(3.21) 



> 



x(2E^(a)-E { J\a)^ + 



2£ 2 (l-b- 1 ) 2 
U-2a(l + e) 



e(£-\ Xl -x 2 \) 



\xi - x 2 \ - £) 



(3.22) 



\x\ - x 2 \ 

To conclude the proof, we need the following lemma. 

Lemma 4 (Calculus lemma general version). Let p 6 i 1 (K+) be non- 
negative, with J p(r)dr = 1. Assume i/iai there are constants a > 0, 6 > 0, c > 0, 
< A < 1 and £ c > swc/i £/ia£ 



(1-A)a + 



f 2 



p{r)dr > I Aa 



p(r)dr 



for all £ > £ c . Then 



1 . . Ac 

-p{r)dr > 

r b + 2c£ r 



(3.23) 



(3.24) 



Proof. Let f(r) = f Q r p(s)ds, and define < < °o by /(£o) = A. In case £q < 
we have f(£ c ) > A and hence 



-p(r)dr > ^ > 



T %. c £ c 

In particular, (|3.24p holds. We can thus assume £q > £ c . For all £ < £q 7 we have 

(1 - X)f(t) < A(l - f(£)) (3.26) 
and hence (|3.23j) implies that 



b f°° 1 

/(£) > c / -p(r)dr 



£ 2 



(3.27) 
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for all l c < £ < £q. Using integration by parts, 

-p(r)dr > ^ + -f'(r)dr = 
r £ r .h r 



■° i 

-sf(r)dr. 



(3.28) 



The first term on the right can be dropped for a lower bound, and in the integrand 
of the second we can use (|3.27|) . We obtain 



1 . . c 
-p(r)ar > - 
r o 



dr 



°°i^ds = c . 



f'(s) 



{mm{s,£ }-£ c )ds. (3.29) 



For a lower bound, we can restrict the integral on the right to s < £q. Using (|3.25|) 
again, this yields 



1 . . c 
—p(r)dr > - 
r b 



/'(*) 



(s — £ c ) ds 



7 Cf(4) - /(4)) - -Ac \ - P (s)ds 



1 



■° 1 



> - A - 24 / -p{r)dr 



o r 



which proves the lemma. 



(3.30) 
□ 



We apply this lemma to (|3.22p , with p being the probability distribution of 
\xi — Xz\ in the state We conclude that, if \I/ satisfies (|3. 211) . then 



This implies (j3~2"j) . 



4- > > 



X(U - 2a(l + e)) 



z ( i-Vi) a +2£ c (U-2a(l + e)) 



(3.31) 
□ 



4. Many Polarons 
We recall that the Hamiltonian Hrf' for N particles is given by (|3.1[i and that 



its ground state and minimum break-up energy were defined in (11 .5[) and (| 1 . 6|) . The 
analogue of Theorem [5] in the TV-particle case is as follows. 

Theorem 3 (Upper bound on the ./V-polaron radius). For any N > 2 and 

e > there is a constant C e (N) > such that for all < 2a(l + e) < U with 
E^\a) < E^\a) and all states * 



f7-2a(l + e) 



(JV) 



C £ (N)(l + U/a) 



E^\a)-E^\a) 



(N) 



(4.1) 



Since (|4.ip holds for all VP, it can be reformulated as an operator inequality. The 



bound is non-trivial only for states \I/ with {^\H\f\^) < E^ 1 '(a), however, which 



exist since E^\a) < E^'(a) by assumption. For approximate ground states, 
satisfying 

{*|J#°I*> < (1 - A^V) + ^ N \a) 
for some A > 0, (|4.ip gives a uniform upper bound on the radius of the multipolaron 
system. The bound depends only on the value of A and does not explode as U — > U c . 



'AN) 
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Proof. We perform a localization similar to that in the two-polaron case, but with 
\xi — 12 1 replaced by the maximum of \x% — Xj\ over all particle pairs Let ifi 

be given as in (|3.4[) - (|3.5p . for some £ > and b > 1. We shall apply (I3.6[) with 



i = max \Xj — a;,- 



(4.2) 



By the IMS localization formula, 



j>0 \ 



i=l fc>0 



j>0 



v3N 



Moreover, for almost every X G 

AT 



2Ei^wi 5 



i=l fc>0 

which can be bounded as in 



fc>0 



eo>4 W) («) 



on the support of In particular, 

»-2 



2{£-£/bf 



(4.3) 



(4.4) 



(4.5) 



For j > 1, we now proceed with the one-particle localization as in the two- 
polaron case, localizing particle i in a ball of radius b J L centered at Ui, for suitably 
chosen L > 0. More precisely, with x given in (|3.10|) . let 



N 

* JiU (x) = *(i)^(t)(w)- M / 2 nx(^(^ 

i=l 



where we denote u = (ut, . . . , ujv). We have 

ll*Vi(*)|| 
and, again using the IMS formula, 



du ( * , u 



H, 



(JV) 



JV^ 2 



6 2 J\L 2 



(4.6) 



(4.7) 



(4.8) 



All particles are now supported in balls Bi of radius b 3 L centered at Uj. Moreover, 
fj U is nonzero only if max^j. \ui — > fr 7-2 '! — 2VL. In particular, if we draw 
open balls of radius R around all the Uj , the resulting set is not connected as long 
as (JV — l)2i? < W~ 2 £ — 2V L. Hence it is possible to split the particles into two 
clusters, in such a way that the distance between any two balls in different clusters 
is at least as big as 



1 



JV — 1 



(b j - 2 £-2b j L) - 2b j L 



N - 1 



(b j ~ 2 £ - 2NVL) 



(4.9) 



We shall choose L < £/(2Nb 2 ) to make this positive. 

The analogue of Lemma [3] that we need now is the following. 
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Lemma 5. Assume that \P is normalized and supported in the set B\ X • • • X Bn , 
where the Bi are balls of radius R. Assume that, for some 1 < n < N — 1, the 
distances dik between balls Bi for i < n and Bk for k > n + 1 are positive. Then, 
for all e > 0, 



A") 



(N-n), 



n N 

E E 

i—l k— n+1 



U -2a{l + e) 



Xk\ 



2a 

dik 



2a(l 



d lk + 4i? 



(4.10) 



Using the path integral representation of the ground state energy, the proof is 
easy; we refer to [5] for details. 

We relabel the particles such that particles 1 < i < n belong to one cluster, and 
rt + l<fc<7Vto the other. The distances dik are all bounded from below by 
They are also bounded above by £b> on the support of With 



Dj := [(N - 1) 



we conclude that 



>E^\a) + (N-l)\\^\\ 



(b>- 2 t - 2NVL) , tb> 
U-2a(l + e) 



(4.11) 



max^fe \Xi - Xk\ 



+ (N - 1) min 

deDj 



2a + 2a(l + e) ^ _ fc _ 
d 



d + AWL J 



b- 2 > 



2^2 



b 2 7T 



2(1 -Ijbf 



1? 



(4.12) 



Here, we have assumed that the minimum over d in the second line is non-negative, 
in which case we can use that n(N — n) > (N — 1) for 1 < n < N — 1. We can now 
argue as in the two-polaron case that, for suitably chosen L, the sum of the terms 
in the second line is positive for £ large enough. In fact, with L = S£/(N — 1) and 



mm 

b- 2 -25A'/(A r -l)< ( i<Ar-l 



1 

"d 



1 



d + 45 



(which is positive for small enough 5) this is the case for all j > 1 if 

- 2 ' b N 



£>L:= 



an e {N-l) 2 V4(l - b- 1 ) 



2bS 2 



(4.13) 



(4.14) 



Recall that t = max^*; \xi — Xk\- We have thus shown the operator inequality 

_2 



H. 



(AT) 



E (N) _ S(JV) 



+ {N-1) 



U ■ 



(a) 2(1 -£/b) 2 
2ail + %(t-£) 



t) 



t 



(4.15) 



for all £ > £ c . The remainder of the proof is as in the two-particle case, applying 
Lemma 2] to the probability distribution of t in a state □ 



5. The Pekar-Tomasevich approximation 

Let us recall, for the reader's convenience, the definition and basic properties of 
the PT model [5]. There is a non-linear differential-integral variational principle 
associated with the polaron problem, which gives the exact ground state energy in 
the limit a — ¥ oo. This variational problem was investigated in detail by Pekar [17) . 



14 



R.L. FRANK, E.H. LIEB, AND R. SEIRINGER 



Pekar and Tomasevich (PT) [T8 j generalized it to the bipolaron, and the extension 
to ./V-polarons obviously follows from [18j . 

The PT functional is the result of a variational calculation and therefore gives an 
upper bound to the ground state energy E^\a). In order to compute H^Hf), 
one takes a * of the form ip <S> $ where ip € L 2 (R 3N ), $e J, and both ip and $ 
are normalized. For a given ip it is easy to compute the optimum $, and one ends 
up with the functional 

N 



j^Jm?" f^J^ 3N \ x i- x i\ 

■dxdy, (5-1) 



<R3 \x-y\ 



where dX — Y\^. =1 dxk, and 



Ptp(x) = 2_] \tp(xi, . . . ,x, . . . ,xn)\ 2 dxi ■ ■ ■ dxi ■ ■ ■ dxN (5.2) 

JR3(«-1) 

with a; at the ?-th position, and dxi meaning that dxi has to be omitted in the 
product dxk- The ground state energy is 

4^(0) = inf {v$PW\ : / \M 2 dX = l) . (5.3) 

I JR3N ) 

The variational argument above gives the upper bound 

4^(a)<4^(a) =4^1(1) o». (5.4) 

(The equality follows by scaling.) For N — 1 this upper bound is due to Pekar; nu- 
merically, one has £W(a) ~ — (0.109)a 2 [15] . Moreover, the minimization problem 
for £ ( 1 )(a) has a unique minimizer (up to translations), see 

The upper bound for N = 1 was widely understood to be asymptotically exact 
for large a. A proof of this was finally achieved by Donsker and Varadhan [4], using 
the functional integral representation and large deviation theory. Later, this fact 
was rederived in [14] by operator methods, and it was shown that the error was no 
worse than a 9 / 5 for large a. The fact that for fixed ratio v = U/a > and for 
N = 2 one has 

lim a- 2 E^\a) = £^ N) (a = 1) (5.5) 



was first noted in [16] ■ This is also valid for arbitrary N. 

We now address the question of whether the infimum 4^ ( a ) ^ s attained, that 
is, whether there is a minimizer. It was proved in [10] that this is the case provided 
' (a) < £~\f\a). This minimum break-up energy is defined as before by 

£$\a) = ^min^ (4"> (a) + (a)) . 

Our next theorem gives an answer in the case of equality, i.e., when £^\a) = 

4 A V). 



BINDING OF POLARONS AND ATOMS AT THRESHOLD 15 

Theorem 4. Let U c > be such that £^^(a) — £^ (a) and assume that there is 
a sequence U n — ¥ U c such that the inftmum £^ (a) is attained. Then the infimum 
£y (a) is attained. 

It is important in our proof that any U c as in Theorem |4] satisfies U c > 2a, as 
shown in [10] . The key input in our proof is once again a bound on the maximal 
distance between the particles. 

Proposition 2. For any v > 1 and N > 2 there is an £q < 00 and a 5q > such 
that, whenever U > 2av and ipu is a minimizer for £^ (a) , then one has 

0(Ei>i-ay| 2 <*o) tf,dX>6o. (5.6) 

r J 

The proof is similar to that of Theorem [3] and is omitted. We only note that, 
since we assume that there is a minimizer, we can take A = 1 in (|3 . 2 1 1) and therefore 
the simple version (|2.5p of the calculus lemma suffices and yields (15.61) . One also 
needs to linearize the functional, as we did in [S] and as we will do later in the proof 
of Theorem 

Proposition 3. Let U n be as in Theorem^ and let tp n be the corresponding min- 
imizers. Then there is a non-zero ip E H 1 (M. 3N ) and a sequence {a n } C M 3 such 
that a subsequence of t/j n (xi — a n , . . . , xn — an) converges weakly in iJ 1 (R 37V ) to -0. 

Proof. Let 77 : [0, 00) — > [0, 1] be a smooth, compactly supported function with 
rj(t) = 1 for t in a neighborhood of zero. We shall prove the following: There is an 
R > 0, a non-zero u € H 1 (M 3N ) and a sequence {a n } C K 3 such that a subsequence 
of Un(xi, . . . ,xn) = v(J2i,j \xi-Xj\ 2 /R 2 )ip n (xi-a n , . . . ,x N -a n ) converges weakly 
in iJ 1 (R 3Ar ) to u. From this, one easily derives the statement of the lemma. 

We first note that ||u n || < ||^ n || = 1 and £i ||V^„|| 2 < ||Vi^„|| a + CR~ 2 
for some C independent of n. Hence the sequence u n is bounded in H 1 (M 3N ), and 
the assertion will follow by an easy extension of [HI Thm. 8.10], provided we can 
show that u n does not converge to zero in measure. 

To prove this, we use the Euler-Lagrange equation satisfied by the functions tp n , 



N 



1 



1 Kj 

2 



with (j> n = 2y/a p^ n * \x\ 1 and A„ = £\j n (a) - ctfjp^ n {x)\x - y\ 1 p^ n (y) dx dy. 
Multiplying the equation by r](J2i j \ x i — Xj\ 2 / R 2 ) 2 ip n {x\, . . . ,xn) and integrating 
we find that 



N 



[ (\v lUn \ 2 - vWn(^K) dx + uY] [ , M " . dX 

JR3N f— ' J R 3N \Xi - Xj\ 



2\ / ..2 



< {\ n + CR- 2 ) / uidX. (5.7) 

Since lim inf n _).oo A n < £ ^ (a) < we can choose R > sufficiently large such 
that liminf„^ 00 (A„ + CR~ 2 ) < 0. Moreover, by our key estimate (15.61) we know 
that, after increasing R if necessary, we have 

lim inf / u 2 dX > . 
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(In order to apply Proposition [2] we use the fact that U c > 2a, see [10].) From this 
we conclude that the only negative term on the left side of (|5.7[) cannot vanish in 
the limit, i.e., 

N 

liminfV / <f> n (xi)uldX>0. (5.8) 
™^°° t~iJ R3N 

We now use the 'pqr theorem' [12j Ex. 2.22] to conclude from (|5.8|) that {u n } 
does not converge to zero in measure. This theorem states, quite generally, that 
for any l<p<(7<r<oo and for any constant C < oo there are constants 
e,S > such that if a function / satisfies ||/|| p < C, \\f\\ r < C and ||/|| q > C _1 , 
then |{|/| > e}| > 6. 

Returning to our concrete situation, we note that, since {«„} is uniformly 
bounded in H 1 (M 3N ), we known from Sobolev inequalities that it is so in L P (R 3N ) 
as well for any 2 < p < 6N/(3N — 2). To apply the 'pqr theorem' we need to show 
that ||w n ||p is uniformly bounded away from zero for some 2 < p < 6N/(3N — 2). 
First, note that, since (y^/tv) is uniformly bounded in ii^R 3 ), 0„ is uniformly 
bounded in L q (M. 3 ) for any 3 < q < oo. Now we choose C > such that r)(t) = 
for t > C and we put 

g n (X) = 9(CR 2 - \*i ~ *j\ 2 ) Yl M*) ■ 

The sequence {<?„} is uniformly bounded in L q (M. 3N ) for any 3 < q < oo, and by 
(15.81) we have 



< 5 < V / <j>n{x 4 )ul dX = [ g n u 2 n dX < 



l5n||g||^||q' 



where q 1 +q' — 1. This shows that (u n ) is uniformly bounded away from zero in 
L 2q '(R 3N ), and we can choose 3 < q < oo in such a way that 2 < 2q' < 6N/(3N-2). 
Hence the 'pqr theorem' implies that u n does not tend to zero in measure, which 
completes the proof of the lemma. □ 

We now turn to the proof of Theorem dJ An important ingredient, which we 
have already used in [5], is that the energy functional can be linearized. In order 
to state this precisely we define for any pair {tp,<p) E H 1 (M 3N ) n i? 1 (M 3 ) (with ijj 
not necessarily normalized) the energy functional 



7V[lM]:=W (\Vii>\ 2 -Va-<f>(xiM 2 ) dX + Uj2 [ dX 
~~[ jr3« ~r!j jr3« f< - x 3 I 



/ \v<p\ 2 dx f 

lOTT J R 3 J R 3 



\i>(X)\ 2 dX, (5.9) 



The crucial observation is that for any ip and <fi one has 

V v [i>A]>Uf ^m^], (5-10) 
with equality if and only if <f) = 2y/a | a; | 1 * p^/\\M\- We are now ready to give the 

Proof of Theorem [^} Because of Proposition [3l after passing to a subsequence and 
a translation we may assume that {ip n } converges weakly in ff^lR 3 ^) to if> ^ 0. 
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Denoting ip n — ip n — ip and (f> n — \x\ 1 * p^ n we claim that 

VuA^nAn] =V Uc [' t P,<t>n]+'Pu c $ n ,<t>-n]+o(l). (5.11) 

Assuming (|5.11j) for the moment, we now finish the proof of Theorem |U Indeed, in 
view of (|5.10[) , we conclude from (|5.11l) that 

4f (a) + 0(1) > ll^ll 2 ^[HVirV] + Unf VtfPlWJ-^n] + 0(1) . 

Since [\\^ n \\^ 1 '4>n] >^u( a ) an< ^ smce ll^nll 2 = 1 ~ IIV'II 2 + °(1)j we learn that 

4f ) (a) + (i)>^ ) [nv.irV] + (i), 

that is, HV'II^V is an optimizer at U = U c . 

It remains to prove (|5.11j) . The analogous assertion separately for the terms 

J2iLi I l^i^l 2 dX and J2i<j I \t^-x -\ ^ 1S an easy consequence of weak conver- 
gence. Moreover, 

\v^ n \ 2 dx= ( \V4> n \ 2 dx ( |VP0I 2 +/ \V4> n \ 2 dx[ |^„(X)| 2 + (1) 



by the weak convergence of ip„ and the fact that J |V0„| 2 dx is uniformly bounded. 
Thus it remains to prove that 

<p n (xi)i> n i>dX = o(l) (5.12) 

3N 

for each i. To prove this, we let R > and split the integral on the left side 
according to whether \X\oo = max{|a;i|, . .., |a;jv|} is bigger or smaller than R. In 
this way we find that 



(xi)ip n ipdX\ 

r 

< sup Halloa fsupll^HII^ OQX^ > R)\\ + U n OQX^ < R)\\ ) . (5.13) 



Here we also used that ||^|| < 1. The fact that {<t> n } is uniformly bounded in 
L°°(M 3 ) follows from the fact that {^/fhp^} is uniformly bounded in iJ 1 (M 3 ). We 
can make the first term on the right side of (|5.13[) as small as we like by choosing 
R large. On the other hand, for any fixed R, the sequence {■0 n #(|X| oo < R)} 
tends to zero in L 2 (R 3N ) by the Rellich-Kondrashov Theorem [TJ Thm. 8.9], and 
therefore we can make the second term on the right side of (|5.13l) as small as 
we like by choosing n large. This completes the proof of (|5. 12|) and therefore of 
Theorem H □ 
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